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1.  Introduction 


In  the  current  theater  of  operations,  warfighters  encounter  threats  that  must  be  detected  at  a 
standoff  distance.  Some  threats  contain  components  whose  permittivity  contrasts  substantially 
with  that  of  the  emplacement;  such  is  the  case  with  many  threats  that  are  buried.  The  reception 
of  a  subsurface  linear  radar  response  from  an  area  whose  surface  is  otherwise  undisturbed 
indicates  the  presence  of  a  threat.  Others  threats  contain  metal  contacts  and  semiconductor 
junctions  whose  nonlinear  electromagnetic  response  contrasts  with  that  of  the  emplacement;  such 
is  the  case  with  radio  frequency  (RF)  electronics.  The  reception  of  a  nonlinear  radar  response 
from  an  area  that  does  not  otherwise  contain  electronics  indicates  the  presence  of  another  class  of 
threat.  The  combined  radar  is  intended  to  detect  both  types  of  threats,  whether  or  not  they  are 
collocated. 

Often,  threats  contain  dielectric,  as  well  as  electronic  components,  hence  they  will  respond  to 
both  linear  and  nonlinear  excitation.  Either  mode  (linear/nonlinear)  will  detect  the  threat.  By 
switching  between  two  radar  modes,  additional  information  about  the  threat  is  received,  and, 
thus,  the  probability  that  it  is  detected  is  improved. 

1.1  Linear  Radar  Versus  Nonlinear  Radar 

Linear  radar  is  well  suited  to  the  detection  of  a  target  whose  complex  permittivity  s  contrasts 
greatly  with  that  of  its  surroundings: 

s  =  s'-j-s",  (1) 

where  s'  is  the  “real”  part,  and  s"  is  the  “imaginary”  part  of  the  permittivity.  The  permittivity 
of  a  material  relative  to  that  of  free  space  is  its  dielectric  constant  sr : 


.  £  t  -  ft 

J  —  =  £r-J-£r 


(2) 


e0  =  8.854  -10“12F/ 

0  /m 


(3) 


Let  the  target  be  illuminated  by  a  radar  wave  as  illustrated  in  figure  1,  and  let  the  electric  field  of 
that  wave  Em  be  represented  by  a  single-tone  sinusoid  of  frequency  fo  and  amplitude  Ed. 

Ein(t)  =  E0COS(27r-f0-t).  (4) 
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Figure  1.  Incident  and  reflected  radar  waves  for 
propagation  normal  to  target. 

Assuming  normal  incidence  (i.e.,  the  direction  of  propagation  of  the  wave  is  normal  to  the 
boundary  of  the  target),  the  reflected  wave  is  (7) 


E,A‘)  = 

\El\cos(27r-f0-t  +  <f>{El}), 

(5) 

—  V-F 

(6) 

—  1  £Lq 

The  frequency  of  the  reflected  wave  is  the  same  as  that  of  the  incident  wave,  but  its  amplitude  is 
scaled  by  the  reflection  coefficient  T  .  For  sr  =  1 ,  T  =  0 ,  the  target  is  transparent  to  the  radar 

wave  traveling  in  air,  and  there  is  no  reflection. 

As  the  contrast  in  sr  between  a  target  and  that  of  the  medium  through  which  the  radar  wave  is 

propagating  increases,  the  strength  of  the  radar  reflection  from  that  target  increases.  The  value 
of  s"  for  a  typical  conductor  (e.g.,  aluminum)  is  greater  than  107  S/m  (2).  Thus,  conductive 

targets  are  very  detectable,  even  if  they  are  buried  or  obscured  by  insulators.  For  insulators,  s"r 
is  near  zero,  but  s'r  can  take  on  a  wide  range  of  values,  from  s'  &  1  for  dry  foam  (3)  up  to 
s'  =  80  for  distilled  water  (4).  Thus,  insulating  targets  are  not  as  detectable,  as  their  radar 
reflections  depend  much  more  strongly  on  s' . 

Nonlinear  radar  exploits  a  completely  different  phenomenon:  it  relies  on  the  nonlinear  properties 
of  a  target  to  convert  a  portion  of  the  incident  radar  wave  into  a  reflected  wave  at  a  different 
frequency.  Most  materials  found  in  nature  are  linear  (with  the  exception  of  ferromagnetics), 
whereas  many  man-made  materials  are  nonlinear  (5-11).  Semiconductor  devices,  such  as  radios 
and  cell  phones,  are  highly  nonlinear. 

A  simple  model  for  RF  nonlinearity  is  the  memoryless  power  series  given  by  (12,  13) 

E,m  (<)  =  (/)  +  a,_El  (/)  +  „,£?(/)  +  ...  =  £ a,E-  (t) ,  (7) 

n= 1 

where  a„  are  complex  power-series  coefficients,  and  Erea  is  the  electric  field  reflected  by  the 
target.  The  value  of  a\  is  the  linear  reflected  response  of  the  target,  T  ;  the  values  \ci2,  a?,, . . . } 
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depend  upon  the  nonlinear  properties  of  the  target.  If  a  nonlinear  target  is  illuminated  by  the 
radar  wave  given  by  equation  4,  the  reflected  wave  is 

oo 

E„«  (0 = Z  K  M2* ■ M  ■  fo  ■'+«>{£„}) .  (8) 

M= 1 

77  _  X  1  ( 2k+M-2)  ® 2k+M-2  j?2k+M-2  /q\ 

k- 1  ^2k+M-3  ^0  ’  \yJ 

k= 1  Z 

which  is  a  sum  of  sinusoids  at  harmonics  M  of fo,  each  with  amplitude  \Em\  and  phase  </>{Em}-  If 
the  radar  measures  EM  =  0  for  all  M>  1,  then  no  nonlinear  target  is  detected.  If  the  radar 
measures  Em  for  some  M>  1,  however,  a  nonlinear  target  is  detected. 

A  combined  radar  detects  targets  using  linear,  as  well  as  nonlinear  reflective  responses.  The 
linear  radar  detects  targets  whose  permittivity  contrasts  with  that  of  the  background,  whereas  the 
nonlinear  radar  detects  targets  whose  electromagnetic  properties  produce  a  change  in  frequency 
between  the  incident  and  reflected  waves. 

1.2  Implementation 

Linear  radar  can  be  implemented  in  different  ways,  which  are  commonly  designated  by  the 
transmit  waveform,  such  as  continuous-wave  (CW),  pulsed  single-tone,  or  chirp.  To  achieve  an 
ultrawide  bandwidth  for  ground  penetration,  as  well  as  an  imaging  resolution,  the  U.S.  Army 
Research  Laboratory  (ARL)  designed  the  Synchronous  Impulse  Reconstruction  (SIRE)  radar 
(14).  The  SIRE  radar  uses  a  single-cycle  impulse  waveform,  two  transmit  antennas,  16  receive 
antennas,  and  multiple  data  traces  collected,  whereas  the  radar  platform  is  in  motion  in  order  to 
form  high-resolution  images  of  surface  and  shallow-buried  targets.  A  single-cycle  impulse  and 
its  spectrum  are  illustrated  in  figure  2a.  An  alternative  design  that  allows  for  more  flexibility  in 
the  transmitted  band  is  the  stepped-frequency  waveform  illustrated  in  figure  2b. 
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Both  impulse  and  stepped- frequency  waveforms  are  broadband.  For  the  impulse,  the  peak  power 
is  high,  but  the  average  power  is  low.  For  the  stepped-frequency  signal,  the  peak  power  and  the 
average  power  are  the  same.  Either  waveform  will  provide  linear  detection  and  ranging. 

One  advantage  of  a  stepped-frequency  design,  however,  is  that  its  underlying  switched- 
frequency  signal  source  is  likely  able  to  dwell  on  a  single  frequency  for  a  long  period  of  time. 

As  dwell  time  increases  while  transmitting  the  same  average  power  in  a  tone  or  a  series  of  tones, 
the  side  lobes  caused  by  interrupting  the  transmission  (e.g.,  turning  the  source  off  or  switching  to 
another  tone)  diminish.  This  extended  dwell  time  is  necessary  in  order  to  minimize  reflected 
linear  side-lobes  from  nonlinear  reflections,  which  are  usually  very  weak. 

Nonlinear  radar  can  also  be  implemented  in  different  ways.  One  popular  technique  is  to  transmit 
a  single  frequency  fo  and  receive  the  target  response  at  the  second  harmonic  of  the  transmitted 
tone,  2fo  (15-19).  A  slight  variation  of  this  technique  tracks  a  Doppler  shift  at  2fo  for  moving 
targets  (20).  Other  variations  chirp  (21)  or  digitally  modulate  (12)  the  transmit  waveform  for 
greater  noise  rejection.  Another  common  technique  is  to  transmit  two  tones  f\  and/2  and  receive 
the  intermodulation  tones  2f\-fi  and  2f2-f\  (6,  22-25).  A  technique  recently  developed  at  ARL 
transmits  at  least  two  tones  and  receives  not  only  a  harmonic  of  the  transmitted  tones  (e.g.,  2f\ 
and  2/2)  but  also  the  mixing  products  of  those  tones  near  that  harmonic  (e.g.,  3/1— /2,  /i  + f 2,  3/2— /i ) 
(26).  Figure  3a  shows  an  example  of  transmit  and  receive  spectra  for  a  nonlinear  radar  that 
transmits  one  tone  and  receives  harmonics  of  that  tone  (27,  28).  Figure  3b  shows  an  example  of 
spectra  for  a  radar  that  transmits  two  tones  and  receives  harmonics,  as  well  as  mixing  products 
near  those  harmonics. 
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Figure  3.  Nonlinear  radar  waveforms:  (a)  single-tone  CW  and  (b)  two-tone  CW. 

In  this  report,  an  experimental  transceiver  is  developed  to  demonstrate  (a)  a  transmit  chain, 
which  generates  waveforms  that  are  appropriate  for  both  linear  and  nonlinear  modes  of 
operation,  and  (b)  a  switchable  receiver  chain,  which  captures  either  linear  or  nonlinear 
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responses  from  a  radar  target.  The  response  to  be  exploited  by  the  nonlinear  receiver  is  the 
second  harmonic  of  the  transmitted  waveform. 

First,  four  different  types  of  waveforms  are  selected  that  are  appropriate  for  linear  and  nonlinear 
radar.  Next,  a  bench-top  architecture  for  the  combined  radar  is  proposed.  Last,  the  experimental 
architecture  is  constructed,  and  measurements  using  a  simulated  (wireline)  radar  channel  are 
performed  to  demonstrate  the  functionality  of  the  linear  and  nonlinear  radar  modes. 


2.  Transmit  Waveforms 


The  four  waveforms  selected  for  the  linear/nonlinear  transmitter  are  the  single-tone  pulse,  the 
multitone  pulse,  the  linear  frequency-modulated  (FM)  chirp,  and  a  stepped- frequency  pulse. 

2.1  Single-Tone  Pulse 

A  mathematical  representation  for  a  single-tone  pulse  produced  by  an  arbitrary  waveform 
generator  (AWG)  is 

L\VG  =  4wC0S(2^-  ./pulse  ^  M0>  (io) 

with  a  carrier  frequency  /pu]SC.  The  amplitude  Aenv  is  computed  from  the  power  of  the  envelope  of 
the  pulse  Pe nv  (in  decibels  referenced  to  1  mW)  by 

Amv  =  ^10«ra/10-2(50  Q)(l0"3  W/mW) .  (11) 

The  pulse  modulation  is  given  by  the  switching  waveform  s(t ): 

s{t)  =  u(t)-u(t-DcT)  =  s(t  +  T)  D  T  -  Tenv ,  (12) 

which  has  a  period  T  and  a  duty  cycle  Dc.  The  pulse  is  active  during  the  time  interval  Tem.  An 
example  of  an  RF  pulse  generated  by  a  Tektronix  AWG7052  is  given  in  figure  4.  A  MATLAB* 
function,  which  generates  a  single-tone  RF  pulse  is  given  in  appendix  A. 

It  should  be  noted  that  (a)  all  signals  in  this  report  were  captured  in  time  by  a  Lecroy 
Wavemaster  8300A  oscilloscope  and  in  frequency  by  an  Agilent  N9342C  spectrum  analyzer;  (b) 
the  sampling  rate  of  the  8300A  oscilloscope  was  20  GS/s,  and  the  resolution  bandwidth  of  the 
N9342C  analyzer  was  1  kHz;  and  (c)  the  amplitude  of  each  waveform  is  less  than  Ae nv  computed 
by  equation  1 1  due  to  the  loss  introduced  by  the  8-ft  RG-58  Subminiature  Version  A  (SMA) 
cable,  which  feeds  each  of  the  signal  capture  instruments. 


MATLAB  is  a  registered  trademark  of  The  MathWorks,  Inc. 
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Figure  4.  Single-tone  RF  pulse  output  by  AWG:  /J,uise  =  900  MHz,  Pcnv  =  0  dBm,  renv  =  1  ps,  Dc  =  10%. 

2.2  Multitone  Pulse 

If,  instead  of  a  single  RF  carrier  frequency,  multiple  frequencies  are  active  during  the  pulse,  a 
multitone  pulse  is  generated: 

VAWG  =  Acne  {C0S  (2^  '  /i‘  t)  +  COS  {lit  ■  f2  •  t)  +  ...  +  COS  {in  •  fN  ■  t)}  J  (t)  ,  (13) 


which  contains  vV  frequencies  given  by/1,  /2,  . .  .  /a?.  In  this  representation,  the  amplitude  of  each 
tone  is  At0 ne,  and  each  tone  begins  at  a  common  initial  phase  (for  maximum  peak-to-average 
ratio,  which  generates  a  maximum  nonlinear  response).  Also,  the  tones  are  centered  at/c  and 
separated  by  /spa ce: 


iV 


/=1 


space 


(14) 


The  active  tones  are  again  modulated  by  the  on/off  pulse  waveform  s(t).  An  example  of  a 
multitone  pulse  is  shown  in  figure  5.  A  MATLAB  function,  which  generates  this  waveform,  is 
given  in  appendix  B. 
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Figure  5.  Multitone  RF  pulse  output  by  AWG:  N=  2  tones,  fe  =  890  MFIz,  PKmc  =  -6  dBm  per  tone, 
renv  =  2  ps,  Dc  =  20%. 

2.3  Linear  Frequency-Modulated  Chirp  Pulse 

A  pulse  whose  carrier  frequency  begins  at  /start  and  increases  linearly  to  /CIK|  over  the  time  interval 
Te nv  is  given  by 

VAWG  =  Anv  C0S  \}K  ■  (/start  +  (^/2)0  '  {]S  (0  k=  (/end  “  /start  )// nv  »  (15) 


where  k  is  the  linear  chirp  rate,  and  Aenv  is  the  amplitude  of  the  pulse  envelope.  An  example  of  a 
linear  FM  chirp  pulse  is  shown  in  figure  6.  A  MATLAB  function,  which  generates  this 
waveform,  is  given  in  appendix  C. 


Frequency  (MHz) 


Figure  6.  Linear  FM  chirp  pulse  output  by  AWG:  /jutrt  =  860  MHz,  /ald  =  900  MHz,  Peav  =  -3  dBm,  Tem  =  4  ps, 
Dc  =  50%. 
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2.4  Stepped-Frequency  Pulse 

A  chirp  whose  carrier  frequency  steps  between  discrete  values  can  be  represented  by 


f{t)  = 


v  =  A 

v  AWG  t 

mvCOs[2 

(16) 

./start 

0  <  t  <  At 

/start  +  A f 

At  <  t  <  2  At 

T 

/start  +2A/ 

2  At  <t<  3  At  At  = 

^env 

N ,  ’ 

(17) 

steps 

.  /end  “A/ 

T  -At  <t<T 

env  env 

where  Aste ps  is  the  number  of  steps,  Tcm  is  the  length  of  the  stepped-frequency  waveform,  Aem  is 
the  amplitude,  A/is  the  spacing  in  frequency  between  each  step,  and  At  is  the  spacing  in  time 
between  each  step.  It  should  be  noted  that  this  representation  for  the  waveform  is  not  phase- 
continuous,  that  is,  the  phase  of  the  waveform  changes  abruptly  across  each  frequency  transition. 


An  example  of  a  stepped-frequency  pulse  is  shown  in  figure  7.  A  MATLAB  function,  which 
generates  this  waveform,  is  given  in  appendix  D. 


Figure  7.  Stepped-frequency  pulse  output  by  AWG:  fstait  =  870  MHz,  /end  =  890  MHz,  A f=  1  MHz, 
P env  =  0  dBm,  renv  =  2.5  ps,  Dc  =  25%. 


3.  Experimental  Transceiver 


An  architecture  having  components  common  to  both  linear  and  nonlinear  modes  for  transmitting 
and  receiving  radar  waveforms  is  necessary  to  minimize  the  size,  weight,  and  power  of  the 
combined  radar  system.  A  bench-top  architecture  for  a  combined  radar  transceiver  is  given  in 
figure  8. 
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Figure  8.  Experimental  combined  radar  architecture. 

In  this  setup,  the  Tektronix  AWG7052  is  the  signal  generator  for  both  linear  and  nonlinear  radar 
waveforms.  The  MiniCircuits  NLP- 1000+  low-pass  filters  are  highly  linear  with  a  passband 
below  1  GHz;  these  filters  remove  much  of  the  transmitter  (Tx)-generated  nonlinear  (harmonic) 
distortion.  The  amplifier  research  AR4W1000  amplifier  boosts  the  power  of  the  AWG  signal  to 
a  level  sufficient  to  excite  nonlinear  responses  from  electronic  targets.  The  HP  778D  dual- 
directional  coupler  provides  one  port  for  sampling  the  forward  (transmit,  Tx)  waveform  and 
another  port  for  sampling  the  reverse  (receive,  Rx)  signal.  The  “Simulated  Radar  Environment” 
consists  of  100  ft  of  SMA  cable  (four  25-ft  cables  in  series),  terminated  by  an  SMA- 
connectorized  target. 

Two  receive  chains  are  selected  by  a  pair  of  Hittite  HMC784MS8GE  switches.  Each  switch  is 
powered  by  5  V  from  the  6-Y/5-A  port  on  an  Agilent  E363 1 A  supply  and  controlled  by 
5  V/0  V  from  the  +25-V/1-A  port.  In  figure  8,  the  “Linear  Rx”  chain  is  selected,  and  the  signal 
is  passed  directly  to  the  8300A  oscilloscope  through  an  SMA  cable.  Alternatively,  the 
“Nonlinear  Rx”  chain  may  be  selected.  Along  the  nonlinear  receiver  path,  the  signal  is  filtered 
by  three  MiniCircuits  VHF-1320+  high-pass  filters  (passbands  above  1.32  GHz,  to  remove  the 
linear  response  from  capture  and  processing)  and  amplified  by  two  MiniCircuits  PSA-5453+  and 
one  MiniCircuits  PSA-545+.  Each  amplifier  is  mounted  on  an  evaluation  board  and  powered  by 
3  V  from  another  E3631A  supply. 

3.1  RF  Signal  Generation  and  Capture 

As  measured  by  an  Agilent  N9923A  network  analyzer  and  observed  in  figure  9,  the  AR4W1000 
provides  more  than  a  40-dB  gain  to  the  transmit  signal.  For  nonlinear  (harmonic)  responses, 
each  NLP- 1000+  filter  attenuates  Tx-generated  distortion  at  frequencies  above  1500  MHz  by 
more  than  40  dB. 
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Figure  9.  Transmitter  amplifier  and  low-pass  filters. 


As  seen  in  figure  10,  the  Tx  and  Rx  coupling  from  the  778D  is  approximately  -20  dB.  Also,  the 
nonlinear  Rx  chain  (measured  from  one  HMC784  “RF  common”  port  to  the  other)  passes  signals 
to  the  8300A  with  a  gain  of  approximately  40  dB,  whereas  the  linear  Rx  chain  passes  signals 
through  with  a  loss  under  3  dB. 


Frequency  (MHz) 


Frequency  (MHz) 


Figure  10.  Directional  coupler  and  linear/nonlinear  receiver  chain. 

3.2  MATLAB  Graphical  User  Interface 

The  AWG7052  generator,  8300A  oscilloscope,  and  E3631A  supplies  are  controlled  via  the 
General  Purpose  Interface  Bus  (GPIB).  Communication  is  established  using  the  Instrument 
Control  Toolbox  in  MATLAB  (v7. 0.0. 19920,  R14).  The  graphical  user  interface  (GUI)  in  figure 
1 1  was  created  using  MATLAB’s  “guide”  function.  The  script  and  functions  that  govern  the 
operation  of  the  GUI  are  given  in  appendices  E  through  K. 
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Figure  11.  GUI  for  experimental  combined-radar  system. 

Using  the  upper  panel,  the  four  different  waveforms  presented  in  section  2  may  be  uploaded  to 
the  AWG  for  transmission  to  the  target.  Using  the  lower  panel,  the  signal  from  the  target  may  be 
captured  using  the  linear  or  nonlinear  receive  chain  and  processed  accordingly. 

For  the  single-tone  pulse,  the  user  may  choose  the  power  of  the  RF  pulse  while  it  is  active  (Pew,), 
the  RF  frequency  (fpuise ),  the  time  interval  during  which  the  pulse  is  active  (Any),  and  the  pulse 
duty  cycle  (Dc).  The  “waveform  name”  is  the  designation  that  appears  on  the  AWG  after  the 
waveform  is  uploaded  and  is  selected  for  waveform  playback. 

For  the  linear  FM  chirp  pulse,  the  user  may  choose  the  power  of  the  chirp  envelope  (Any),  the 
frequency  at  which  the  chirp  starts  (fs tart),  the  frequency  at  which  the  chirp  ends  (/cnd),  the  time 
interval  during  which  the  frequency  linearly  changes  from /sta rt  to  /enci  (renv),  and  the  duty  cycle  of 
the  waveform  (Dc). 

For  the  multitone  pulse,  the  user  may  choose  the  number  of  tones  (N),  power  per  tone  (Aone),  the 
frequency  at  which  the  tones  are  centered  (fc ),  the  time  interval  during  which  the  pulse  is  active 
(Tenv),  and  the  pulse  duty  cycle  (Z)c).  The  frequency  separation  between  the  tones  (/space)  is 
automatically  set  to  \ITsm,  so  that  the  shortest  waveform  necessary  to  achieve  N,fc,  and  Te nv  with 
negligible  frequency  aliasing  is  uploaded  to  the  AWG. 

For  the  stepped- frequency  waveform,  the  user  may  choose  the  power  of  the  pulse  envelope 
(Pem),  the  frequency  at  which  the  stepping  starts  (/tart),  the  frequency  at  which  the  stepping  ends 
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(/end)?  the  time  interval  during  which  the  frequency  steps  from^tart  to  fen a  (renv),  the  step  size  (A/), 
and  the  duty  cycle  of  the  waveform  ( Dc ). 

In  agreement  with  figure  8,  the  sampled  Tx  signal  is  fed  to  channel  2  of  the  8300A  oscilloscope, 
and  the  sampled  Rx  signal  is  fed  (through  the  linear/nonlinear  receive  chain)  to  channel  3.  The 
user  chooses  the  voltage  scale  per  channel,  the  total  data  collection  time  per  trace,  and  the 
number  of  integrations  (i.e.,  the  number  of  data  traces  averaged  before  capture). 

The  user  chooses  the  trigger  level  and  source.  In  the  experimental  setup,  Marker  1  from  the 
AWG7052  is  fed  to  the  External  trigger  port  on  the  8300A. 

The  user  chooses  the  receiver  (Rx)  mode  and  types  a  name  for  the  native  MATLAB  (MAT)  file 
that  will  store  the  time-sampled  Tx  and  Rx  voltage  vectors. 

Upon  pressing  the  “Upload. . .”  button  inside  of  one  of  the  upper  subpanels,  the  appropriate 
waveform  is  generated  and  sent  to  the  AWG.  A  new  figure  panel  (not  shown)  appears,  which 
plots  the  software-generated  waveform  in  frequency  and  time  to  confirm  that  the  signal  the  user 
intended  has  been  uploaded. 

Upon  pressing  the  “Capture. . .”  button  inside  the  lower  subpanel,  the  corresponding  signal 
received  from  the  target  is  recorded  by  the  oscilloscope  and  processed  in  MATLAB.  A  second 
figure  panel  (shown  in  section  4)  appears,  which  plots  the  raw  Tx  and  Rx  data  in  time.  A  third 
figure  panel  (also  shown  in  section  4)  appears,  which  plots  the  correlation  of  the  Tx  and  Rx 
voltage  samples. 


4.  Wireline  Experiments 


Several  experiments  were  conducted  in  order  to  demonstrate  the  performance  of  the  bench-top 
combined-radar  transceiver:  three  different  waveforms  (pulse,  linear  chirp,  stepped-frequency), 
two  Rx  modes  (linear  and  nonlinear),  two  linear  targets  (open-circuit,  matched  load),  and  one 
nonlinear  target  (Family  Radio  Service  [FRS]  radio). 

4.1  Linear  Rx,  Chirp  Waveform,  Open-Circuit  Target 

Figure  12  shows  the  result  of  the  linear  data  capture  and  processing  when  reflecting  a  chirp  from 
an  open-circuit  target.  Figure  12a  plots  the  raw  Tx  and  Rx  data.  Figure  12b  plots  the  cross 
correlation  of  the  Tx  and  Rx  signals: 

[ym  *  V,« ] W  =  J"  I (<) ■ •  ('  +  O' * ,  (18) 

where  time  has  been  mapped  to  distance  using  the  velocity  of  propagation  of  an  RF  signal  in  the 
MiniCircuits  CBL  25-ft  coaxial  lines  (dielectric  constant  sv  «  2.1): 
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A  factor  of  1/2  is  used  in  equation  19,  because  the  distance  plotted  is  half  the  round-trip  distance 
from  the  transmitter  (i.e.,  from  the  coupler  output  port)  to  the  target  (i.e.,  to  the  end  of  the  100-ft 
coaxial  line)  to  the  receiver  (i.e.,  back  to  the  coupler  output  port). 
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Figure  12.  Radar  data,  chirp  Tx  waveform,  linear  Rx  mode,  open-circuit  target:  fst art  =  880  MHz,/cnd  =  920  MHz, 
Pcnv  =  0  dBm,  Tcnv  =  1  ps,  Dc  =  10%,  (a)  raw  Tx  and  Rx  data  and  (b)  correlation  of  Tx  and  Rx  waveforms. 


Cross  correlation  is  a  basic  form  of  target  ranging.  The  peak  of  Ftrans*  Viec  (as  a  function  of 
distance)  indicates  the  distance  from  the  transmitter  to  the  target. 

In  figure  12a,  a  relatively  constant-amplitude  pulse  is  visible  in  the  sampled  Tx  channel,  and  a 
distorted  pulse  is  visible  in  the  sampled  Rx  channel.  These  waveforms  are  expected,  given  the 
frequency-dependent  characteristic  of  the  coupler  in  figure  10.  In  figure  12b,  a  sine  function  is 
visible  along  with  several  sidelobes.  This  shape  is  expected  from  the  cross  correlation  of  two 
chirps.  The  peak  of  the  sine  function  is  visible  at  a  distance  of  d  =  101  ft.  The  calculated 
distance  is  very  close  to  the  length  of  the  coaxial  line  (and  slightly  higher  because  the  calculation 
does  not  account  for  the  length  of  the  Rx  chain). 

4.2  Linear  Rx,  Chirp  Waveform,  Open  Circuit  Versus  Matched  Load 

Figure  13a  shows  the  result  of  the  cross  correlation  when  reflecting  a  chirp  with  a  wider 
bandwidth  than  that  of  section  4. 1  from  an  open-circuit  target.  Figure  13b  is  the  correlation 
when  reflecting  the  same  waveform  from  a  matched  (50-Q)  load.  Two  results  are  notable:  (1) 
the  peak  is  sharper  when  the  bandwidth  of  the  Tx  waveform  is  wider  and  (2)  very  little  signal 
reflects  from  the  matched  load.  Both  results  are  expected  and  indicate  proper  operation  of  the 
transmitter  and  the  linear  receive  chain. 
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(a) 


(b) 


Figure  13.  Radar  data,  chirp  Tx  waveform,  linear  Rx  mode:  /start  =  860  MHz,  /end  =  940  MHz,  Pe nv  =  0  dBm, 
renv  =  1  ps,  Dc  =  10%,  (a)  open-circuit  target,  (b)  matched-load  target. 


4.3  Nonlinear  Rx,  Pulse  Waveform,  Nonlinear  Target 

Figure  14  shows  the  result  of  the  nonlinear  data  capture  when  reflecting  an  RF  pulse  from  a 
nonlinear  target:  a  Motorola  T4500  radio  whose  antenna  has  been  replaced  by  an  SMA  end- 
launch  connector.  A  13-dB  attenuator  is  placed  between  the  end  of  the  coaxial  line  and  the  FRS 
radio.  Figure  14a  plots  the  raw  Tx  and  Rx  data  along  a  5-ps  time  scale.  Figure  14b  plots  the 
same  raw  data  along  a  2-ns  time  scale  between  t  =  500  ns  and  t  =  502  ns. 


Figure  14.  Radar  data,  RF  pulse  Tx  waveform,  linear  Rx  mode,  FRS  radio  target:  /puise  =  900  MHz,  Penv  =  0  dBm, 
Tcnv  =  1  ps,  Dc  =  10%,  (a)  raw  data,  complete  time  scale  and  (b)  raw  data,  zoomed-in  time  scale. 
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It  is  not  obvious  from  figure  14a  that  the  transceiver  is  detecting  the  nonlinear  response  from  the 
target.  In  figure  14b,  however,  the  observed  response  is  clearly  nonlinear,  because  the  frequency 
of  the  received  signal  (1800  MHz)  is  twice  that  of  the  transmitted  signal  (900  MHz). 

4.4  Nonlinear  Rx,  Stepped-Frequency  Waveform,  Nonlinear  Target 

Nonlinearity  is  also  visible  in  the  frequency  domain  when  the  Tx  and  Rx  signals  are  captured 
with  a  spectrum  analyzer.  Figure  15  provides  such  captures  for  a  stepped-frequency  waveform 
and  the  FRS  radio  target.  The  signal  output  from  the  AWG  (and  filtered  by  a  NLP- 1000+)  is 
plotted  above,  and  the  received  spectrum  is  plotted  below.  For  Fawg,  all  of  the  spectral  content 
is  centered  at /=  900  MHz,  and  no  spectral  content  exists  near  2/=  1800  MHz.  For  Pie c,  all  of 
the  spectral  content  is  centered  at  2/=  1 800  MHz,  and  no  spectral  content  exists  near /=  900 
MHz. 


Figure  15.  Tx  and  Rx  frequency  content,  stepped-frequency  Tx  waveform,  nonlinear  Rx  mode, 

FRS  radio  target:  fsXmt  =  890  MHz,  /cnd  =  910  MHz,  A f=  1  MHz,  Pcnv  =  0  dBm, 

Tcnv  =  2  ps,  De  =  20%. 

4.5  Nonlinear  Rx,  Chirp  Waveform,  Nonlinear  Target 

Figure  16  shows  the  result  of  the  nonlinear  data  capture  and  processing  when  reflecting  a  chirp 
from  the  FRS  radio.  Figure  16a  plots  the  raw  Tx  and  Rx  data.  Figure  16b  plots  the  cross 
correlation  of  the  Tx  and  Rx  signals: 

[L  *  rj(0  =  J“  VL.  (<) ' •  V,K  (t  +  0 dr ,  (20) 

where  the  Tx  signal  used  for  the  correlation  is  a  filtered  second  harmonic  of  the  captured  Limns: 
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and  /?bpf  is  a  bandpass  fdter  with  passband  edges  /(  =  3/c/2  and  /u  =  5/c/2,  where/  =  (/start  + 

/end)/2. 

A  sine  function  is  again  visible  and  centered  at  =  103  ft.  This  distance  is  longer  than  d  =  101  ft 
measured  previously,  because  the  nonlinear  Rx  chain  contains  slightly  more  propagation  delay 
(through  the  filters  and  amplifiers)  than  the  linear  Rx  chain  (SMA  cable,  pass  through). 
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Figure  16.  Radar  data,  linear  chirp  Tx  waveform,  nonlinear  Rx  mode,  FRS  radio  target:  /tart  =  880  MHz, 
d  =  920  MHz,  Pem  =  0  dBm,  rCTV  =  1  ps,  Dc  =  10%,  (a)  raw  Tx  and  Rx  data  and  (b)  correlation  of  Tx  and  Rx 
waveforms. 


4.6  Nonlinear  Rx,  Chirp  Waveform,  Nonlinear  Versus  Open-Circuit  Target 

The  nonlinear  Rx  chain  was  tested  against  a  purely  linear  target  in  order  to  demonstrate  that  the 
transceiver  does  not  indicate  detection  if  the  target  is  linear,  and  the  radar  is  listening  for  a 
nonlinear  response.  Figure  17  gives  the  result  of  this  test,  which  is  performed  with  a  chirp 
waveform. 

From  figure  17a,  it  is  clear  that  the  radar  registers  a  detection  (d  =  103  ft)  when  the  target  is 
nonlinear,  and  the  Rx  is  expecting  a  nonlinear  response.  From  figure  17b,  it  is  clear  that  the 
nonlinear  Rx  chain  does  not  register  a  detection  when  the  target  is  linear. 
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Figure  17.  Radar  data,  chirp  Tx  waveform,  nonlinear  Rx  mode:  /star t  =  890  MHz,^nd  =  910  MHz,  PenY  =  0  dBm, 
renv  =  1  jlxs,  Dc  =  10%,  (a)  FRS  radio  target  and  (b)  open-circuit  target. 


5.  Conclusions 


A  combined-radar  transceiver  has  been  constructed,  which  enables  basic  target  ranging  in  both 
linear  and  nonlinear  (harmonic)  receive  modes.  The  transceiver  was  constructed  using  an 
arbitrary  waveform  generator  as  the  signal  source,  a  high-speed  digitizing  oscilloscope  as  the 
signal  capture  device,  and  commercial  off-the-shelf  (COTS)  components  for  the  radar  front-end 
(amplification,  filtering,  and  switching).  A  100-ft  length  of  SMA  cable  terminated  in  an  open 
circuit  simulated  a  linear  radar  target;  the  same  cable  terminated  in  an  SMA-connectorized  FRS 
radio  simulated  a  nonlinear  radar  target.  A  MATLAB  GUI  was  developed  in  order  to  control  the 
transceiver  remotely.  The  associated  script  and  helper  functions  are  provided  in  the  appendices. 
Ranging  to  the  target  was  demonstrated  experimentally  using  RF  pulses,  linear  FM  chirps,  and 
stepped-frequency  waveforms. 
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Appendix  A.  MATLAB  Function  for  Generating  Single-Tone  RF  Pulsest 


t  This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  [v,  f_sample]  =  awg_pulse (  f_pulse,  T_env,  P_env,  duty_cycle,  f_sample_max 


f  res  =  1/T  env; 

%  frequency  resolution  (Hz) 

f  sample  =  2*f  res; 

while  (f  sample  <=  f  sample  max/2) 
f  sample  =  2*f  sampler- 

end 

%  adjust  AWG  sample  rate 

f  pulse  =  round(f  pulse/f  res)*f  res; 

%  adjust  pulse  carrier  frequency 

t  sbb  =  1/f  sample; 
t  =  0 : t  sbb : T  env-t  sbb; 

A  =  sqrt(10A(P  env/10)  *2*50*10^-3)  ; 

%  sampling  time 
%  time  vector 

%  amplitude  of  pulse  when  ON 

x  =  A  *  cos  (2*pi*f  pulse. *t) ; 

%  sinusoid  equation 

if  (duty  cycle  <=  0)  | |  (duty  cycle  >=  1) 

v  =  x; 

else 

%  continuous  wave 

v  =  [x  zeros (1,  (round (l/duty_cycle) -1) *length (x) )]  ;  %  pulsed 

end 
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Appendix  B.  MATLAB  Function  for  Generating  Multitone  RF  Pulsest 


t  This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  [v,  f  sample]  =  awg  multitone (  N, 

f  c,  T  env,  P  tone,  duty  cycler  f  sample  max  ) 

f  res  =  1/T  env; 
f  space  =  2*f  res; 

%  resolution  based  on  envelope  length 
%  tone  spacing 

f  sample  =  2*f  res; 

while  (f  sample  <=  f  sample  max/2) 
f  sample  =  2*f  sampler- 

end 

%  adjust  AWG  sample  rate 

f  c  =  round(f  c/f  res)*f  res; 

%  adjust  center  frequency 

A  =  sqrt (10A (P_tone/10) *2*50*10^-3)  ; 

%  amplitude  of  each  tone 

t  sbb  =  1/f  sample; 
t  =  0 : t  sbb : T  env-t  sbb; 

%  sampling  time 
%  time  vector 

xi  =  zeros  (1, length (t) ) ; 
xq  =  zeros (1, length (t) ) ; 

%  initialize  baseband  I/Q  vectors 

%  generate  multitones  with  in-phase  and  quadrature  components 
for  n  =  1 : N/2 

xi  =  xi  +  (A) *cos (2*pi*f_space* (abs (n) -1/2) *t)  +  (A) *cos (2*pi*f_space* (abs (n) -1/2 ) *t ) 
xq  =  xq  +  (A)  *sin  (2*pi*f_space*  (abs  (n)  -1/2)  *t)  -  (A)  *sin  (2*pi*f_space*  (abs  (n) -1/2  )  *t ) 

end; 

x  =  xi . *cos (2*pi*f_c*t  )  -  xq. *sin (2*pi*f_c*t  ); 


if  (duty  cycle  <=  0)  | |  (duty  cycle  >=  1) 

v  =  x; 

else 

%  continuous  wave 

v  =  [x  zeros (1, (round (l/duty_cycle) -1) ^length (x) )] ;  %  pulsed 

end 
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Appendix  C.  MATLAB  Function  for  Generating  RF  Chirp  Pulses§ 


§ 


This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  [v,  f_sample]  =  awg_chirp (  f_start,  f_end,  T_env,  P_env,  duty_cycle,  f_sample_max 


f  res  =  1/T  env; 

%  frequency  resolution  (Hz) 

f  sample  =  2*f  res; 

while  (f  sample  <=  f  sample  max/2) 
f  sample  =  2*f  sampler- 

end 

%  adjust  AWG  sample  rate 

f  start  =  round(f  start/f  res)*f  res; 
f  end  =  round (f  end/ f  res)*f  res; 

%  adjust  start  frequency 
%  adjust  end  frequency 

t  sbb  =  1/f  sample; 
t  =  0 : t  sbb : T  env-t  sbb; 

A  =  sqrt(10/v(P  env/10)  *2*50*10^-3)  ; 

%  sampling  time 
%  time  vector 

%  amplitude  of  chirp  when  ON 

k  =  (f  end-f  start) /T  env; 

%  chirp  rate  (Hz/s) 

x  =  A  *  cos(2*pi*(f  start! (k/2 ) *t) . *t) ; 

%  chirp  equation 

if  (duty  cycle  <=  0)  | |  (duty  cycle  >=  1) 

v  =  x; 

else 

%  continuous  wave 

v  =  [x  zeros (1, (round (l/duty_cycle) -1) ^length (x) )] ;  %  pulsed 

end 
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Appendix  D.  MATLAB  Function  for  Generating  Stepped-frequency  RF 

Pulses** 


This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  [v,  f_sample]  =  awg_stepped(  f_start ,  f_end, f_delta, T_env, P_env, duty_cycle, 
f_res  =  1/T_env;  %  frequency  resolution  (Hz) 

f_sample  =  2*f_res; 

while  (f_sample  <=  f_sample_max/2 )  %  adjust  AWG  sample  rate 

f_sample  =  2*f_sample; 

end 


freq  =  ( f_start : f_delta : f_end-f_delta) ;  %  instantanous  frequencies 

N  —  length ( freq) ; 


%  adjust  instantaneous  frequencies 
for  counter  =  1 : length ( freq) 

freq (counter )  =  floor (freq (counter) /f_res) *f_res; 

end 


t_sbb  =  l/f_sample; 
t  =  0 : t_sbb : T_env-t_sbb; 

A  =  sqrt (10A (P_env/10) *2*50*10^-3)  ; 


%  sampling  time 
%  time  vector 

%  amplitude  of  chirp  when  ON 


index  =  ceil  (  (1 : length (t) )/ (length (t) /N) ) ; 
f  =  freq(index); 
x  =  A  *  cos (2*pi* (f . *t) ) ; 


if  (duty_cycle  <=  0)  | |  (duty_cycle  >=  1) 

v  =  x;  %  continuous  wave 

else 

v  =  [x  zeros (1,  (round (l/duty_cycle) -1) ^length  (x) )] ;  %  pulsed 

end 


sample_max  ) 
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Appendix  E.  MATLAB  Script  for  the  Graphical  User  Interfaced 


tt 


This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  varargout  =  gremlin (varargin) 


gui_Singleton  =  1; 
gui_State  =  struct (' gui_Name ' , 

' gui_Singleton ' , 

' gui_OpeningFcn ' 

' gui_OutputFcn ' , 

' gui_LayoutFcn ' , 

' gui_Callback ' , 
if  nargin  &&  ischar (varargin { 1 } ) 

gui_State . gui_Callback  =  str2func (varargin { 1 }) ; 

end 


mfilename,  . . . 
gui_Singleton,  . . . 
@gremlin_OpeningFcn, 
@gremlin_OutputFcn, 

[] )  ; 


if  nargout 

[varargout { 1 : nargout } ]  =  gui_mainf cn (gui_State ,  varargin {:}) ; 

else 

gui_mainfcn (gui_State,  varargin { : } ) ; 

end 


%%%  executes  before  the  GUI  is  made  visible  %%% 


function  gremlin_OpeningFcn (hOb j ect ,  eventdata,  handles,  varargin) 

%  connect  to  the  DC  power  supplies  (Agilent  E3631A) 
global  dca  deb; 

dca  =  visa ( ’ ni ' , 'GPIBO: :5: : INSTR ' ) ; 
deb  =  visa ( 'ni ', 'GPIBO :: 4 :: INSTR' ) ; 

f open (dca) ; 
fopen (deb) ; 


fprintf (dca, 

[ ' APPLy 

P6V, 

5  V'] ) ; 

% 

apply 

5V 

to 

RF  switches 

fprintf (deb, 

[ 'APPLy 

P6V, 

3  V']); 

% 

apply 

3V 

to 

low-noise 

amplifiers 

fprintf (dca, 

[ 'APPLy 

P25V, 

0  V«l); 

% 

start 

up 

in 

linear  Rx 

mode 

fprintf (deb, 

[ 'APPLy 

P25V, 

5  V' ] ) ; 

fprintf (dca,  [ ' OUTPut : STATe  ON ' ] ) ; 
fprintf (deb,  [' OUTPut : STATe  ON ' ] ) ; 
f close (dca) ; 
fclose (deb) ; 

%  connect  to  the  AWG  (Tektronix  AWG7052) 

%  and  to  the  oscilloscope  (Lecroy  Wavemaster  8300A) 
global  awg  osc; 

awg  =  visa ( 'ni 1 , 'GPIBO : :2 :: INSTR' ) ; 
osc  =  visa ( 'ni ',  'GPIBO :: 6 :: INSTR' ) ; 

fopen (awg) 

fwrite (awg,  ' AWGCONTROL : D0UTPUT1 : STATE  1');  %  turn  'direct  output'  on 

fclose (awg) 


%  define  strings  for  the  Tx 
%  Tx  channel  averaged  (F2) 
global  data_trc_l  data_trc_ 
data_trc_l  =  ' C2 ' ; 
data_trc_2  =  ' C3 '  ; 
avg_trc_l  =  ' F2 '  ; 
avg_trc_2  =  ' F3 ' ; 
sweeps  =  1; 


channel  (C2),  Rx  channel  (C3) , 
and  Rx  channel  averaged  (F3) 
avg_trc_l  avg_trc_2  sweeps; 


%  turn  C2  and  C3  off;  turn  F2  and  F3  on 
fopen (osc) 

fprintf (osc, [avg_trc_l  ' : DEF  EQN, ' ' AVG ( '  data_trc_l  '), SUMMED '', SWEEPS , '  num2str (sweeps )] ) 
fprintf (osc,  [avg_trc_2  '  : DEF  EQN,  ''AVG ('  data_trc_2  '), SUMMED '', SWEEPS ,  '  num2str ( sweeps )] ) 
fprintf  (osc,  [data_trc_l  '  : TRACE  OFF' ] ) ; 
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fprintf (osc, [data_trc_2 
fprintf (osc, [avg_trc_l  ' 
fprintf (osc, [avg_trc_2  ' 
fprintf (osc, [avg_trc_l  ' 
fprintf (osc, [avg_trc_2  ' 
f close (osc) 


: TRACE  OFF ' ] ) ; 
TRACE  ON ' ] ) ; 
TRACE  ON ' ] ) ; 
FRST ' ] ) ; 

FRST ' ] ) ; 


reset  the  averaged  sweeps 


handles . output  =  hObject; 
guidata (hObj ect ,  handles); 


generate  single-tone  RF  pulse 


function  pushbuttonl_Callback (hOb j ect ,  eventdata,  handles) 
global  awg; 

%  user-defined  waveform  name  to  be  uploaded  to  the  AWG 
global  wf_name; 

wf_name  =  get (handles . edit 6 ,  'String'); 


%  user-defined  pulse  power  while  the  pulse  is  active 
global  P_env; 

P_env  =  str2num (get (handles . editl ,  'String')); 


%  user-defined  carrier  frequency  for  the  RF  pulse 
global  f_pulse; 

f_pulse  =  str2num (get (handles . edit2 , 'String')); 


%  user-defined  time  for  the  RF  pulse  to  be  active 
global  T_env; 

T_env  =  str2num (get (handles . edit4 , 'String')); 


%  user-defined  duty  cycle 
global  duty_cycle; 

duty_cycle  =  str2num (get (handles . edit5 ,' String '))/100; 


%  user-defined  maximum  sampling  rate  for  the  AWG 
global  f_sample_max; 

f_sample_max  =  str2num (get (handles . edit3 ,' String ')) ; 

%  generate  the  single-tone  pulse  waveform 
%  and  compute  an  appropriate  AWG  sampling  frequency 
global  v  fsample; 

[v,f_sample]  =  awg_pulse (f_pulse, T_env, P_env, duty_cycle, f_sample_max) ; 
awg_plot (v, f_sample) 


global  samples; 
samples  =  length (v) ; 


%  provide  two  markers  (triggers)  in  sync  with  the  AWG  output 

marker_l  =  zeros (1, samples ) ; 

marker_l (1 : round (samples/10) )  =  1; 

marker_2  =  zeros (1, samples ) ; 

marker_2 (1 : round (samples/10) )  =  1; 


global  buffer; 
buffer  =  samples; 
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set  (awg,  ' InputBuf f erSize ' , buf f er ) ; 
set (awg,  ' OutputBuff erSize ' , buf f er )  ; 


%  convert  the  real  waveform  data  to  binary 
%  and  upload  the  binary  data  to  the  AWG  for  playback 
[binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2); 
awg_upload (awg,  buffer,  wf_name,  binblock,  header,  bytes); 


set  the  AWG  sampling  frequency, 
and  turn  the  AWG  output  on 


fopen (awg) ; 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
f close (awg) ; 


' AWGC : D0UT1  ON'); 

[ ' S0URCE1 : FREQUENCY 
[ ' SOURCE 1 : WAVEFORM  " 
'OUTPUT  ON'  )  ; 

' AWGCONTROL : RUN ' )  ; 


select  the  waveform  to  be  played, 


num2str  (f_ 
wf  name  '  ' 


sample/lCU 

'  ]); 


9) 


GHZ '  ]  )  , 


global  f_trans; 
f_trans  =  f_pulse; 


%%%  generate  linear  RF  chirp  pulse  %%% 


function  pushbutton2_Callback (hOb j ect ,  eventdata,  handles) 
global  awg; 

%  user-defined  waveform  name  to  be  uploaded  to  the  AWG 
global  wf_name; 

wf_name  =  get (handles . editl2  ,  'String'); 


%  user-defined  chirp  power  while  the  pulse  is  active 
global  P_env; 

P_env  =  str2num (get (handles . edit8 , 'String')); 


%  user-defined  start  frequency  for  the  chirp 
global  f_start; 

f_start  =  str2num (get (handles . edit 9 ,' String ')) ; 

%  user-defined  end  frequency  for  the  chirp 
global  f_end; 

f_end  =  str2num (get (handles . editl3 ,' String ')) ; 

%  user-defined  length  of  the  chirp 
global  T_env; 

T_env  =  str2num (get (handles . editlO ,' String '))  ; 

%  user-defined  duty  cycle 
global  duty_cycle; 

duty_cycle  =  str2num (get (handles . editll ,' String ' ) ) / 1 0  0 ; 


%  user-defined  maximum  sampling  rate  for  the  AWG 
global  f_sample_max; 

f_sample_max  =  str2num (get (handles . edit3 , 'String')); 


%  generate  the  RF  chirp  waveform 

%  and  compute  an  appropriate  AWG  sampling  frequency 
global  v  f_sample; 

[v, f_sample]  =  awg_chirp (f_start, f_end, T_env, P_env, duty_cycle, f_sample_max) 
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awg_plot (v, f_sample) 


global  samples; 
samples  =  length (v) ; 


%  provide  two  markers  (triggers)  in  sync  with  the  AWG  output 

marker_l  =  zeros (1, samples ) ; 

marker_l (1 : round (samples/10) )  =  1; 

marker_2  =  zeros (1, samples) ; 

marker_2 (1 : round (samples/10) )  =  1; 


global  buffer; 
buffer  =  samples; 

set  (awg,  ' InputBuf f erSize ' ,buffer) ; 
set (awg, ' OutputBuff erSize ' ,buffer) ; 


%  convert  the  real  waveform  data  to  binary 
%  and  upload  the  binary  data  to  the  AWG  for  playback 
[binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2); 
awg_upload (awg,  buffer,  wf_name,  binblock,  header,  bytes); 


set  the  AWG  sampling  frequency, 
and  turn  the  AWG  output  on 


fopen (awg) ; 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
f close (awg) ; 


' AWGC : D0UT1  ON  '  )  ; 

[ ' S0URCE1 : FREQUENCY 
[ ' SOURCE 1 : WAVEFORM  " 
'OUTPUT  ON'  )  ; 

' AWGCONTROL : RUN ' )  ; 


select  the  waveform  to  be  played. 


num2str  (f_ 
wf  name  '  ' 


sample/lCU 

']); 


9) 


GHZ '  ]  )  , 


global  f_trans; 

f_trans  =  (f_start  +  f_end)/2; 


%%%  generate  multitone  RF  pulse  %%% 


function  pushbutton4_Callback (hOb j ect ,  eventdata,  handles) 
global  awg; 

%  user-defined  waveform  name  to  be  uploaded  to  the  AWG 
global  wf_name; 

wf_name  =  get (handles . editl8 ,' String ') ; 
global  N; 

N  =  str2num (get (handles . editl 9 ,  'String')); 


%  user-defined  power  per  tone  while  the  pulse  is  active 
global  P_tone; 

P_tone  =  str2num (get (handles . editl4 ,' String ')) ; 

%  user-defined  center  frequency  for  the  tones 
global  f_c; 

f_c  =  str2num (get (handles . editl5 , 'String')); 


%  user-defined  time  for  the  multitone  pulse  to  be  active 
global  T_env; 

T_env  =  str2num (get (handles . editl 6 ,' String ')) ; 
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%  user-defined  duty  cycle 
global  duty_cycle; 

duty_cycle  m  str2num (get (handles . editl7 ,  ' String ')) /100 ; 

%  user-defined  maximum  sampling  rate  for  the  AWG 
global  f_sample_max; 

f_sample_max  =  str2num (get (handles . edit3 ,  ' String ')) ; 

%  generate  the  multitone  waveform 

%  and  compute  an  appropriate  AWG  sampling  frequency 
global  v  fsample; 

[v, f_sample]  =  awg_multitone (N, f_c, T_env, P_tone, duty_cycle, f_sample_max) 
awg_plot (v, f_sample) 


global  samples; 
samples  =  length (v); 


%  provide  two  markers  (triggers)  in  sync  with  the  AWG  output 

marker_l  =  zeros (1, samples ) ; 

marker_l (1 : round (samples/10) )  =  1; 

marker_2  =  zeros (1, samples) ; 

marker_2 (1 : round (samples/10) )  =  1; 


global  buffer; 
buffer  =  samples; 

set (awg, ' InputBuf f erSize ' ,buffer) ; 
set (awg,  ' OutputBuff erSize ' ,buffer)  ; 


%  convert  the  real  waveform  data  to  binary 
%  and  upload  the  binary  data  to  the  AWG  for  playback 
[binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2); 
awg_upload (awg,  buffer,  wf_name,  binblock,  header,  bytes); 


set  the  AWG  sampling  frequency, 
and  turn  the  AWG  output  on 


fopen (awg) ; 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
f close (awg)  ; 


' AWGC : D0UT1  ON'); 

[ ' S0URCE1 : FREQUENCY 
[ ' SOURCE 1 : WAVEFORM  " 
'OUTPUT  ON' ) ; 

' AWGCONTROL : RUN ' )  ; 


select  the  waveform  to  be  played. 


num2str  (f_ 
wf  name  '  ' 


sample/10/ 

'  ]  )  ; 


9) 


GHZ ' ] ) , 


global  f_trans; 
f_trans  =  f  c; 


%%%  generate  stepped-f requency  pulse  %%% 


function  pushbutton5_Callback (hOb j ect ,  eventdata,  handles) 
global  awg; 

%  user-defined  waveform  name  to  be  uploaded  to  the  AWG 
global  wf_name; 

wf_name  =  get (handles . edit2 4 ,' String ') ; 
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%  user-defined  power  at  each  frequency  while  the  pulse  is  active 
global  P_env; 

P_env  =  str2num (get (handles . edit2 0 ,  'String')); 


%  user-defined  start  frequency  for  the  steps 
global  f_start; 

f_start  =  str2num (get (handles . edit2 6 ,  ' String ')) ; 

%  user-defined  end  frequency  for  the  steps 
global  f_end; 

f_end  =  str2num (get (handles . edit2 1 ,  ' String ')) ; 

%  user-defined  frequency  step  size 
global  f_delta; 

f_delta  =  str2num (get (handles . edit25 ,  ' String ')) ; 

%  user-defined  time  for  the  stepped  waveform  to  be  active 
global  T_env; 

T_env  =  str2num (get (handles . edit22 ,  ' String '))  ; 

%  user-defined  duty  cycle 
global  duty_cycle; 

duty_cycle  =  str2num (get (handles . edit23 ,  ' String ')) /100 ; 

%  user-defined  maximum  sampling  rate  for  the  AWG 
global  f_sample_max; 

f_sample_max  =  str2num (get (handles . edit3 ,  ' String ')) ; 

%  generate  the  stepped-f requency  waveform 
%  and  compute  an  appropriate  AWG  sampling  frequency 
global  v  f_sample; 

[v, f_sample]  =  awg_stepped (f_start, f_end, f_delta, T_env, P_env, duty_cycle, f_sample_max) 
awg_plot (v, f_sample) 


global  samples; 
samples  =  length (v); 


%  provide  two  markers  (triggers)  in  sync  with  the  AWG  output 

marker_l  =  zeros (l,samples) ; 

marker_l (1 : round (samples/10) )  =  1; 

marker_2  =  zeros (l,samples) ; 

marker_2 (1 : round (samples/10) )  =  1; 


global  buffer; 
buffer  =  samples; 

set (awg, ' InputBuf f erSize ' , buf f er ) ; 
set (awg,  ' OutputBuff erSize ' , buf fer )  ; 


%  convert  the  real  waveform  data  to  binary 
%  and  upload  the  binary  data  to  the  AWG  for  playback 
[binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2); 
awg_upload (awg,  buffer,  wf_name,  binblock,  header,  bytes); 


set  the  AWG  sampling  frequency, 
and  turn  the  AWG  output  on 


fopen (awg) ; 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 
fwrite (awg, 


' AWGC : D0UT1  ON ' ) ; 

[ ' S0URCE1 : FREQUENCY 
[ ' SOURCE 1 : WAVEFORM  " 
'OUTPUT  ON' ) ; 

' AWGCONTROL : RUN ' ) ; 


select  the  waveform  to  be  played. 


num2str ( f_sample/10A 9 )  '  GHZ']); 

wf_name  '"']); 
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f close (awg) ; 


global  f_trans; 

f_trans  =  (f  start  +  f_end)/2; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  clear  waveforms  &  reset  AWG  %%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


function  pushbutton3_Callback (hOb j ect ,  eventdata,  handles) 

global  awg; 

f open (awg) ; 
fwrite (awg, ' *RST ' ) ; 
fwrite (awg, ' *CLS ' ) ; 
pause  ( 1 ) 

fwrite (awg,  ' AWGCONTROL : D0UTPUT1 : STATE  1'); 
f close (awg) ; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  capture  and  process  Rx  waveform  %%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


function  pushbutton7_Callback (hOb j ect ,  eventdata,  handles) 


global  mode; 

switch  get (handles . popupmenu3 , 'Value') 
case  1 

mode  =  1;  %  select  LINEAR  receive  mode 

case  2 

mode  =2;  %  select  NONLINEAR  receive  mode 

otherwise 

mode  =  1 ; 

end 

global  dca  deb; 
fopen (dca) ; 
f open (deb) ; 
if  (mode  ==  2) 


fprintf (dca, 

[ ' APPLy 

P25V, 

5 

V'  ]  )  ; 

% 

set  RF  switches  to 

fprintf (deb, 

else 

[ 'APPLy 

P25V, 

0 

V’]  )  ; 

% 

nonlinear  Rx  chain 

fprintf (dca, 

[ 'APPLy 

P25V, 

0 

V  ]); 

% 

set  RF  switches  to 

fprintf (deb, 

end 

[ 'APPLy 

P25V, 

5 

V  ]); 

% 

linear  Rx  chain 

f close (dca) ; 
f close (deb) ; 

global  awg  osc; 

global  trig_src; 

switch  get (handles . popupmenu2 , 'Value') 
case  1 


trig  sre  =  'EX'; 
case  2 

trigger  =  external 

trig  sre  =  ' Cl ' ; 
case  3 

% 

trigger  =  channel  Cl 

trig  sre  =  ' C4 ' ; 
otherwise 

trig  sre  =  ' EX ' ; 

4 

trigger  =  channel  C4 

end 
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%  user-defined  trigger  level 
global  trig_lev; 

trig_lev  =  str2num (get (handles . edit37 ,  ' String ')) *10A-3 ; 

%  user-defined  total  time  for  Tx  and  Rx  data  to  be  collected 
global  t_total; 

t_total  =  str2num (get (handles . edit32 ,  ' String ')) *10A-6; 
t_div  =  t  total/10; 


%  user-defined  voltage  per  division  visible  on  the  oscilloscope 
global  v_div_l  v_div_2; 

v_div_l  =  str2num (get (handles . edit33 ,  ' String ')) *10A-3  ; 
v_div_2  m  str2num  (get  (handles  .  edit34  ,  '  String  '))  *10A-3 ; 


%  user-defined  voltage  offset  on  the  oscilloscope 
global  v_offset; 
v_offset  =  Oe-3; 


%  set  trigger  delay  to  1/2  of  the  data  collection  time  so  that 
%  the  maximum  amount  of  data  is  collected  AFTER  the  trigger  point 
global  trig_delay; 
trig_delay  =  -t_total/2; 


%  user-defined  Matlab  file  to  store  Tx  and  Rx  time  &  voltage  vectors 
global  data_name; 

data_name  =  get (handles . edit3 6 ,  ' String ' ) ; 


global  trig_mode; 

trig_mode  =  'SINGLE’;  %  use  a  single  trigger  for  each  sweep 

global  f_capture; 
f_capture  =  20e9; 


global  buffer; 

buffer  =  10*t_div*f_capture  +  26; 
set (osc, ' InputBuf f erSize ' , buf f er ) ; 
set (osc, ' OutputBuff erSize ' ,buffer) ; 

global  data_trc_l  data_trc_2  avg_trc_l  avg_trc_2  sweeps; 
sweeps  =  str2num (get (handles . edit38 ,  'String')); 


%  set  the  oscilloscope  to  record  data  according  to  user-input  values 
fopen (osc) ; 

[ ' TRSE  EDGE, SR, '  trig_src  ' , HT, OFF ' ] ) ; 

[ ' TRMD  '  trig_mode] ) ; 

[  '  TRDL  '  num2str  (trig_delay/10/"-6)  '  US  ]); 

[ ' EX : TRLV  '  num2str ( trig_lev/10 A-3 )  '  MV' ] ) ; 

[ ' TDIV  '  num2str(t  div/lCT-6)  '  US ' ] ) ; 


fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fclose (osc) ; 


[data_ 

[data_ 

[data_ 

[data_ 


trc 

trc 

trc 

trc 


VDIV  '  num2str (v_div_l/10A-3)  '  MV' ] ) ; 

VDIV  '  num2str (v_div_2/10 A-3 )  '  MV']); 

OFST  '  num2str (v_off set/10 A-3 )  '  MV']), 

OFST  '  num2str (v  off set/10 A-3 )  '  MV']), 


%  turn  the  appropriate  data  traces  on/off 

%  and  set  two  of  the  traces  to  each  average  the  Tx  or  Rx  data  stream 
fopen (osc) ; 

fprintf (osc, [data_trc_l  ': TRACE  OFF']); 
fprintf  (osc,  [data_trc_2  ': TRACE  OFF ' ] ) ; 
fprintf (osc, [avg_trc_l  ': TRACE  ON']); 
fprintf (osc, [avg_trc_2  ': TRACE  ON']); 
fprintf (osc, [avg_trc_l  ' : FRST ' ] ) ; 
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fprintf (osc, [avg_trc_2  ' : FRST '  ]  )  ; 

fprintf (osc, [avg_trc_l  '  :  DEF  EQN, 1 ' AVG ( '  data_trc_l  '), SUMMED '', SWEEPS ,  '  num2str (sweeps) ] ) 
fprintf (osc, [avg_trc_2  ' : DEF  EQN, ' ' AVG ( '  data_trc_2  '), SUMMED '', SWEEPS , '  num2str (sweeps) ] ) 
for  counter  =  1 : sweeps 
fprintf (osc, ' *TRG ' ) ; 
fprintf (osc, 'WAIT' ) ; 

end 

f close (osc) ; 


%  download  the  oscilloscope  data  into  Matlab 

[v_trans,  v_rec,  t_trans,  t_rec]  =  osc_capture (  osc,  buffer,  avg_trc_l,  avg_trc_2  ); 


%  stop  the  waveform  generator  playback 
f open (awg) ; 

fwrite (awg,  'OUTPUT  OFF'); 
f write (awg,  ' AWGCONTROL : STOP ' )  ; 
f close (awg) ; 


global  f_trans; 


%  process  and  display  the  Tx  and  Rx  data 

osc_process (  v_trans,  v_rec,  t_trans,  t_rec,  mode,  f_trans  ); 

%  save  the  raw  Tx  and  Rx  vectors  to  a  MAT  file 
save (data_name, ' t_trans ' , ' v_trans ' , ' t_rec ' , ' v_rec ' ) ; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%  initialize  all  text  boxes,  buttons,  etc.  %%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  varargout  =  gremlin_OutputFcn (hOb j ect ,  eventdata,  handles) 
varargout{l}  =  handles . output ; 

function  editl_Callback (hOb j ect ,  eventdata,  handles) 

function  editl_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) ; 

end 

function  edit2_Callback (hOb j ect ,  eventdata,  handles) 

function  edit2_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) ; 

end 

function  edit3_Callback (hOb j ect ,  eventdata,  handles) 

function  edit3_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) ; 

end 


function  edit4_Callback (hOb j ect ,  eventdata,  handles) 
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function  edit4_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set  (hOb j  ect,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, ’BackgroundColor’ ,  get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit5_Callback (hOb j ect ,  eventdata,  handles) 

function  edit5_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, ’BackgroundColor’ , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit 6_Callback (hObj ect ,  eventdata,  handles) 

function  edit 6_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ’ ,  ' white ' ) ; 

else 

set (hObject, 'BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit8_Callback (hOb j ect ,  eventdata,  handles) 

function  edit8_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit 9_Callback (hObj ect ,  eventdata,  handles) 

function  edit 9_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ’ ,  ' white ' ) ; 

else 

set (hObject, 'BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editlO_Callback (hObj ect ,  eventdata,  handles) 

function  editlO_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, 'BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editll_Callback (hObj ect ,  eventdata,  handles) 

function  editll_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, 'BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editl2_Callback (hOb j ect ,  eventdata,  handles) 
function  editl2_CreateFcn (hOb j ect ,  eventdata,  handles) 
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if  ispc 

set  (hObj ect,  1 BackgroundColor '  ,  ' white ' ) ; 

else 

set  (hOb j ect,  ' BackgroundColor ' , get ( 0 ,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editl3_Callback (hOb j ect ,  eventdata,  handles) 

function  editl3_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editl4_Callback (hObj ect ,  eventdata,  handles) 

function  editl4_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editl5_Callback (hOb j ect ,  eventdata,  handles) 

function  editl5_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editl 6_Callback (hOb j ect ,  eventdata,  handles) 

function  editl 6_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editl7_Callback (hObj ect ,  eventdata,  handles) 

function  editl7_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hOb j  ect ,  ' BackgroundColor ' , get ( 0 ,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editl8_Callback (hObj ect ,  eventdata,  handles) 

function  editl8_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  editl 9_Callback (hObj ect ,  eventdata,  handles) 
function  editl 9_CreateFcn (hObj ect ,  eventdata,  handles) 
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if  ispc 

set  (hObj ect,  1 BackgroundColor '  ,  ' white ' ) ; 

else 

set  (hOb j ect,  ' BackgroundColor ' , get ( 0 ,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit2 0_Callback (hOb j ect ,  eventdata,  handles) 

function  edit2 0_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit2 l_Callback (hOb j ect ,  eventdata,  handles) 

function  edit2 l_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit22_Callback (hOb j ect ,  eventdata,  handles) 

function  edit22_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit23_Callback (hObj ect ,  eventdata,  handles) 

function  edit23_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit2 4_Callback (hOb j ect ,  eventdata,  handles) 
function  edit2 4_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit25_Callback (hObj ect ,  eventdata,  handles) 

function  edit25_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, 'BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit2 6_Callback (hOb j ect ,  eventdata,  handles) 
function  edit2 6_CreateFcn (hOb j ect ,  eventdata,  handles) 
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if  ispc 

set  (hOb j ect,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, ’BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit32_Callback (hOb j ect ,  eventdata,  handles) 

function  edit32_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hObj ect,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  ’BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit33_Callback (hOb j ect ,  eventdata,  handles) 

function  edit33_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit34_Callback (hObj ect ,  eventdata,  handles) 

function  edit34_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit3 6_Callback (hOb j ect ,  eventdata,  handles) 

function  edit3 6_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, 'BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit37_Callback (hOb j ect ,  eventdata,  handles) 
function  edit37_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject,  'BackgroundColor' , get (0,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  popupmenu2_Callback (hOb j ect ,  eventdata,  handles) 
function  popupmenu2_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, 'BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 
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function  popupmenu3_Callback (hOb j ect ,  eventdata,  handles) 

function  popupmenu3_CreateFcn (hOb j ect ,  eventdata,  handles) 
if  ispc 

set (hObj  ect,  1 BackgroundColor ' ,  ' white ' ) ; 

else 

set  (hOb j  ect,  ' BackgroundColor ' , get ( 0 ,  ' def aultUicontrolBackgroundColor ' ) ) 

end 

function  edit38_Callback (hOb j ect ,  eventdata,  handles) 

function  edit38_CreateFcn (hObj ect ,  eventdata,  handles) 
if  ispc 

set (hOb j  ect ,  ' BackgroundColor ' ,  ' white ' ) ; 

else 

set (hObject, ’BackgroundColor' , get (0, ' def aultUicontrolBackgroundColor ' ) ) 

end 
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Intentionally  Left  Blank. 
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Appendix  F.  MATLAB  Function  for  an  Ideal  Bandpass  Filter** 


it 


This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  y_out  =  bandpass (  x_in,  y_in,  start_freq,  end_freq  ) 

%  y_out  =  bandpass (  x_in,  y_in,  start_freq,  end_freq  ) 

%  --  applies  an  ideal  bandpass  filter  to  ' y_in '  as  a  function  of  ' x_in' 

%  --  start_freq  =  lower  bandpass  corner  frequency 

%  --  end_freq  =  upper  bandpass  corner  frequency 
%  --  y_out  =  filtered  result,  still  a  function  of  ' x_in ' 

T  =  x_in (length (x_in) )  -  x_in(l);  %  full  period  of  x 

tau  =  x_in(2)  -  x_in(l);  %  delta (x)  sample  length 

Y_temp  =  fft (y_in) ;  %  Fast  Fourier  Transform 

if  mod (length (Y_temp) , 2 )  ==  1 

frequency  =  [  (0 : floor (length (Y_temp) /2 ))* (1/T)  fliplr ( (1 : floor (length (Y_temp) /2) )* (-1/T) ) 

]  ; 

else 

frequency  =  [  ( 0 : length (Y_temp) /2 )*( 1/T)  fliplr ( (1 : length (Y_temp) /2-1) * (-1/T) )  ]; 

end 


for  n  =  1 : length ( frequency ) 

if  (  frequency (n)  >=  -l*start_f req  )  &&  (  frequency (n)  <=  start_freq  )  ... 
|  |  (  frequency (n)  <=  -l*end_freq  )  I  I  (  frequency (n)  >=  end_freq  ) 
Y_temp(n)  =  0; 

end 

end 

y_out  =  real (if ft (Y_temp) ) ;  %  inverse  FFT 
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Appendix  G.  MATLAB  Function  for  Converting  Real  AWG  Waveforms  to 

Binary  Data§§ 


§§  This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  [binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2) 

y  =  2  *  single (v) ; 

samples  =  length (v); 

marker_l  =  uint8 (marker_l) ; 
marker_2  =  uint8 (marker_2 ) ; 


marker_l  =  bitshift (uint8 ( logical (marker_l )), 6 ) ; 
marker_2  —  bitshift (uint8 ( logical (marker_2 )), 7 ) ; 
m  =  marker_l  +  marker_2; 

binblock  =  zeros ( 1 , samples*5 ,  ' uint8 ') ; 
for  k=l: samples 

binblock ( (k-1) *5+1 : (k-1) *5+5)  =  [ typecast (y (k) ,  ' uint8 1 )  m(k)]; 

end 

bytes  =  num2str (length (binblock) ) ; 

header  =  ['#'  num2str (length (bytes ) )  bytes]; 
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Appendix  H.  MATLAB  Function  for  Uploading  Waveform  Data  to  the 

AWG*** 


This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  awg_upload (  awg,  buffer,  wf_name,  binblock,  header,  bytes) 
fopen (awg) ;  %  connect  to  the  waveform  generator 

fwrite (awg, [ ' WLISt : WAVeform: DELete  " '  wf_name  '"']); 

f write (awg, [ 'WLI ST : WAVEFORM: NEW  "'  wf_name  ' " , '  num2str (buf f er )  ' , REAL ' ] ) 

cmd  =  [ 'WLIST : WAVEFORM: DATA  wf_name  header  binblock]; 

bytes  =  length (cmd); 
if  buffer  >=  bytes 
fwrite (awg, cmd) 

else 

awg.EOIMode  =  'off'; 

for  i  =  1 : buf f er : bytes-buf f er 

fwrite (awg, cmd (i : i+buf fer-1 ) ) 

end 

awg.EOIMode  =  'on'; 
i  =  i  +  buffer; 
fwrite (awg, cmd (i : end) ) 

end 

fclose (awg) ; 
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Appendix  I.  MATLAB  Function  for  Plotting  Waveforms  Uploaded  to  the 

AWGttt 


ttt 


This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  awg_plot (  v,  f_sample  ) 


t  =  (0 : length (v) -1) *l/f_sample; 

T  =  t  (length (t) )  +  l/f_sample; 
f  =  -f_sample/2 : 1/T : ( f_sample/2 ) -1/T; 


%  time  vector 
%  total  time  per  waveform 
%  frequency  grid 


V  =  f ftshift (f ft (2* (v) ) /length (v) ) ; 

PdB  =  10*logl0 (abs (V)  . A2/  (2*50*10A-3) ) ; 


%  calculate  FFT 
%  calculate  power  spectrum 


figure ( 1 ) 


%  plot  in  time  &  frequency 


subplot (1,2,1) 
plot (t/10A-6,v/10A-3) 

axis ([-Inf  Inf  1 . l*min (v/10 A-3 )  1 . l*max (v/10 A-3 ) ] ) 

ylabel ( 'Voltage  (mV) ' ) 

xlabel ( ' Time  ( \mus ) ' ) 

subplot (1,2,2) 

plot (f/10A6, PdB) 

axis([0  Inf  -Inf  max (PdB) +10 ] ) 

ylabel ( ' Power  (dBm) ' ) 

xlabel (' Frequency  (MHz)') 
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Appendix  J.  MATLAB  Function  for  Recording  Real  Waveforms  Using  the 

Oscilloscope#* 


***  This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  [v_l,  v_2 ,  t_l ,  t_2 ]  =  osc_capture (  osc,  buffer,  data_trc_l,  data_trc_2  ) 

fopen(osc);  %  connect  to  the  oscilloscope 

fprintf (osc,  'ARM') 

pause ( 1 ) 

fprintf (osc,  [data_trc_l  ' :WF?  ALL']);  %  grab  raw  data  trace 

raw_data  =  f read (osc, buff er,  ' uint8 ') ; 

%  calculate  voltage  amplitude  per  step 

v_quant_temp  =  raw_data ( 181 ) *1 6A 6  +  raw_data ( 180 ) *1 6 A4  +  raw_data ( 17 9 ) *1 6 A2  +  raw_data ( 178 ) *1 6 A 
v_quant  =  typecast (uint32 (v_quant_temp) ,' single ') ; 

%  calculate  voltage  offset 

v_of f set_temp  =  raw_data ( 185 ) *1 6 A 6  +  raw_data (184 ) *16A4  +  raw_data ( 183 ) *1 6 A2  + 
raw_data (182) *1 6 A0 ; 

v_offset  =  typecast (uint32 (v_offset_temp) ,' single ') ; 

%  calculate  time  interval  per  step 

t_quant_temp  =  raw_data (201) *16A6  +  raw_data (200) *16A4  +  raw_data (199) *16A2  +  raw_data ( 1 98 ) *1 6 A 
t_quant  =  typecast (uint32 ( t_quant_temp) , ' single ' ) ; 

%  calculate  time  offset 

t_of f set_temp  =  raw_data (209) *16A14  +  raw_data (208) *16A12  +  raw_data (207 ) *16A10  +  ... 
raw_data (206) *16A8  raw_data (205) *16A6  +  raw_data (204 ) *16A4  +  raw_data (203) *16A2  +  ... 
raw_data (202) *16A0; 

fprintf (osc,  [data_trc_l  ' :WF?  DAT1 ' ] ) ;  %  grab  voltage  trace 

raw_data_v  m  f read (osc, buff er, ' int8 ') ; 

%  convert  raw  data  to  time  &  voltage  waveforms 

v_l  =  v_quant*raw_data_v (23 : length (raw_data_v) -1 )  +  v_offset; 

t_l  =  t_quant* ( 0 : length (v_l ) -1 )  +  t_offset; 

pause ( 1 ) 

fprintf (osc,  [data_trc_2  ' :WF?  ALL']);  %  grab  raw  data  trace 

raw_data  =  f read (osc, buff er,  ' uint8 ') ; 

%  calculate  voltage  amplitude  per  step 

v_quant_temp  =  raw_data ( 181 ) *1 6 A 6  +  raw_data ( 180 ) *1 6 A4  +  raw_data ( 17 9 ) *1 6 A2  +  raw_data ( 178 ) *1 6 A 
v_quant  =  typecast (uint32 (v_quant_temp) ,' single ') ; 

%  calculate  voltage  offset 

v_of f set_temp  =  raw_data ( 185 ) *1 6 A 6  +  raw_data ( 184 ) *1 6 A4  +  raw_data ( 183 ) *1 6 A2  + 
raw_data (182) *1 6 A0 ; 

v_offset  =  typecast (uint32 (v_offset_temp) ,' single ') ; 

%  calculate  time  interval  per  step 

t_quant_temp  =  raw_data (201) *16A6  +  raw_data (200) *16A4  +  raw_data (199) *16A2  +  raw_data ( 1 98 ) *1 6 A 
t_quant  =  typecast (uint32 ( t_quant_temp) , ' single ' ) ; 

%  calculate  time  offset 

t_of f set_temp  =  raw_data (209) *16A14  +  raw_data (208) *16A12  +  raw_data (207 ) *16A10  +  ... 
raw_data (206) *16A8  raw_data (205) *16A6  +  raw_data (204 ) *16A4  +  raw_data (203) *16A2  +  ... 
raw_data (202) *16A0; 

t_offset  =  typecast (uint 64 (t_of f set_temp) ,' double ') ; 

fprintf (osc,  [data_trc_2  ' :WF?  DAT1 ' ] ) ;  %  grab  voltage  trace 

raw_data_v  =  f read (osc, buff er, ' int8 ') ; 

%  convert  raw  data  to  time  &  voltage  waveforms 

v_2  =  v_quant*raw_data_v (23 : length (raw_data_v) -1 )  +  v_offset; 

t_2  =  t_quant* ( 0 : length (v_2 ) -1 )  +  t_offset; 

f close (osc) ; 


0 

0 

0 

0 
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Appendix  K.  MATLAB  Function  for  Processing  Radar  Data  from  the 

Oscilloscope§§§ 


1  This  appendix  is  presented  in  its  original  form  without  editorial  change. 
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function  osc_process (  v_l,  v_2,  t_l,  t_2 ,  mode,  f_trans  ) 

delta_t  =  t_2(2)  -  t_2(l); 

T  =  t_l (length (t_l) )  -  t_l (1) ; 

%  for  nonlinear  mode,  correlate  Rx  with  2nd  harmonic  of  Tx 
if  (  mode  ==  2  ) 

v _ lx  =  v _ 1 . A2 ; 

v_lx  =  bandpass (t_l,v_lx, 1 . 5*f_trans, 2 . 5*f_trans) ; 

else 

v_lx  =  v_l ; 

end 

%  compute  correlation  and  associated  time  vector 
v_corr  =  xcorr ( v_lx, v_2 ) ; 

t_corr  =  delta_t* ( 0 : length (v_corr) -1 )  -  T; 

%  adjust  the  velocity  of  propagation  by  dielectric  constant 
dielectric_PTFE  =  2.1; 

velocity_PTFE  =  0.9836e9  /  sqrt (dielectric_PTFE ) ; 

%  compute  a  distance  vector  based  on  2X  transit  time 
d_corr  =  -1/2  *  t_corr  *  velocity_PTFE ; 

%  find  the  peak  of  the  correlation  waveform 
[v_max, index]  =  max (abs (v_corr) ) ; 
v_max  =  abs (v_corr (index) ) ; 
d_max  =  d_corr (index); 


figure (2)  %  plot  the  raw  Tx  and  Rx  data 

subplot (2,1,1) 
plot  (t_l/10A-6,  v_l,  '  -  '  ) 
axis ( [-Inf  Inf  -1.5  1.5]) 
ylabel ( 1 Tx  Voltage,  Ch2  (V)') 

subplot (2,1,2) 

plot  (t_2/10/"-6,  v_2/10/"-3,  '-') 
axis ([-Inf  Inf  -250  250]) 
ylabel ('Rx  Voltage,  Ch  3  (mV) ') 
xlabel('Time  (\mus) ' ) 

figure (3 ) 
subplot (2,1,1) 

plot (t _ i/10A  — 9,v _ 1,  ) 

axis ( [500  502  -1.5  1.5]  ) 
ylabel ('Tx  Voltage,  Ch2  (V) ') 

grid 

subplot (2,1,2) 

plot (t_2/10A-9,v_2/10A-3,  '-') 
axis ( [500  502  -150  150]  ) 
ylabel ('Rx  Voltage,  Ch  3  (mV) ') 
xlabel('Time  (ns)') 
grid 

figure (5)  %  plot  the  Tx  and  Rx  correlation 

plot (d_corr , v_corr ) 

axis ([50  400  r1000  1000]) 

%  text (250 , 800 ,[' peak  @  {\itd}  =  '  num2str (round (d_max) )  '  ft']); 

ylabel ('{ \bf Correlation,  { \itV}_{ trans }  *  { \itV}_{ rec } } ' ) 
xlabel  ('{ \bf Distance  to  Target}  (ft)') 


%  plot  Tx  and  Rx  data  samples 
%  to  confirm  Rx  frequency 
%  is  twice  the  Tx  frequency 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ARL 

U.S.  Army  Research  Laboratory 

AWG 

arbitrary  waveform  generator 

COTS 

commercial  off-the-shelf 

CW 

continuous-wave 

FM 

frequency  modulated 

FRS 

Family  Radio  Service 

GPIB 

General  Purpose  Interface  Bus 

GUI 

graphical  user  interface 

RF 

radio  frequency 

Rx 

receiver 

SIRE 

Synchronous  Impulse  Reconstruction 

SMA 

Subminiature  Version  A 

Tx 

transmitter 
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